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We explore the centrality dependence of the properties of the dense hadronic matter created 
in y/SNN = 200 GeV Au-Au collisions at RHIC. Using the statistical hadronization model, we fit 
particle yields known for 11 centrality bins. We present the resulting model parameters, rapidity 
yields of physical quantities, and the physical properties of bulk matter at hadronization as function 
of centrality. We discuss the production of strangeness and entropy. 
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I. INTRODUCTION 

The measurement of hadron rapidity yields at the Rel- 
ativistic Heavy Ion Collider (RHIC) facilitates a study of 
the physical properties of the hadronic fireball at time of 
hadronization (i.e., when these particles are produced). 
The objective of this work is to understand the impact 
parameter (reaction volume) dependence of the fireball 
bulk properties. We search for a change of the reaction 
mechanism as function of centrality: if a new state of 
matter is formed in central and semi-central nuclear AA 
collisions, but not in pp reactions, one would naively ex- 
pect a visible change in some physical bulk properties for 
a sufficiently small number of reaction participants. 

We consider, at the top RHIC energy snn = 
200 GeV Au-Au, the 11 centrality bins in which the 
7T^~ , K ± , p and p rapidity yields have been recently pre- 
sented, see table I and table VIII in Ref. 0. These pre- 
cise experimental results involving a full range of central- 
ity motivate this effort. We wish to establish, at the high 
level of precision now available for the RHIC -^/snn = 200 
GeV run, what a rapid change of the particle ratios such 
as K + /tt + , K~ /tt~ as function of centrality means both 
for the bulk physical properties of the fireball, and for 
the statistical hadronization model (SHM) parameters 
dependence on centrality. 

These six particle rapidity yield results are comple- 
mented with STAR results for the ratios K*(892)/K" 0, 
and <f>/K- |j. Both K*(892)/K" and (f>/K~ do not show 
a large centrality dependence, but we make an effort to 
account for any dependence in our analysis. 

We considered the difference between STAR 0, and 
PHENIX (submitted for publication, Q) (^-results. We 
illustrate the situation in Fig.2]||. The lines show our 
best fit results to STAR (top panel), PHENIX (middle 
panel) and combined data set (bottom panel) . The inte- 
grated yields agree for the top two panels with those re- 
ported by the experimental collaborations. We note that 



the integrated yield derived from the combined data fit 
(bottom panel of Fig.QJ, to all available 10% centrality 
^-yields, is not compatible with the PHENIX yield. This 
is so, since the evaluation of the integrated PHENIX <j>- 
yield depends on the lowest m± measured yield. This 
data point appears to be a 1.5 s.d. low anomaly com- 
pared to the many STAR ^-results available at low m±. 
This possibly statistical fluctuation materially influences 
the total integrated PHENIX 0-yield. 

The ratio K*(892)/K~ anchors and confirms the chem- 
ical freeze-out temperature T, which is the only param- 
eter on which this ratio depends. The ratio <^>/K - com- 
prises a multi-strange particle and anchors and confirms 
the chemical conditions at freeze-out. For this reason 
use of these yield results is of essence to obtain the pre- 
cision results we present here. However, in principle 
the resonance yields maybe significantly altered by post- 
hadronization processes 0,0,13) or their observableyield 
could be impacted by decay product rescattering 

The chemical non-equilibrium hadronization model de- 
scribes the experimental data analyzed very well. More- 
over, in this model, one expects near coincidence of the 
thermal and chemical freeze-out. In this limit, there is 
no post-hadronization resonance yield evolution, or sig- 
nificant decay product rescattering. We have never come 
across the need for SHM adjustments of yields of the K*- 
resonance or <f>. Therefore, we do not pursue the develop- 
ment of kinetic yield evolution models for these particles. 
The interested reader can follow up these developments 
in Refs. 



II. STATISTICAL HADRONIZATION MODELS 
AND STATISTICAL PARAMETERS 

The statistical hadronization model is, by definition, a 
model of particle production in which the birth process 
of each particle fully saturates (maximizes) the quantum 
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FIG. 1: (color online) Measured <f>-m± distributions 
dNi/dm±dy for 10% most central collisions at ^/sjvjv = 200 
GeV, lines are best fits. From top to bottom: STAR (sub- 
script S), PHENIX (subscript P) and combined all data (sub- 
script A). 



mechanical probability amplitude, and thus, the yields 
are determined by the appropriate integrals of the acces- 
sible phase space jljj- F° r a system subject to global 
dynamical evolution, such as collective flow, this is un- 
derstood to apply within each local co-moving frame el- 
ement. The results presented here were obtained using 
the numerical package SHARE (Statistical Hadronization 
with REsonances) |T1|. 

The question, if SHM is indeed consistent with the 
wealth of RHIC data available today, comes to mind. 
Our comprehensive study of central reactions at central 
rapidity, for both ^/snn = 130 and 200 GeV, suggests 
so strongly [lj ■ Systematic study of particle production 
for a wide reaction energy range confirms applicability 
of the SHM, for review see ^jjj]. At RHIC there are two 
2 s.d. (standard deviation) exceptions among the model 
agreement with hadron particle yields at RHIC-200: 
1) the (preliminary) H/fi = 1.01 ± 0.08 yield ratio pf. 
with the central value greater than unity while, on gen- 



eral grounds, at finite baryon density this ratio should be 
smaller than unity; 

2) the (preliminary) A++/p = 0.24 ± 0.06 0, Q E3, 
which in statistical hadronization models is half as large. 
We note that the central value of this result means that, 
after removal of descendants from weak decays, nearly all 
protons observed should have been a primary A. 

In addition to the 8 particle (relative) yields consid- 
ered, we also enforce three supplemental constraints: 

1) strangeness conservation, i.e., the (grand canonical) 
count of s quarks in all hadrons equals such s count for 
each rapidity unit; 

2) the electrical charge to net baryon ratio in the final 
state is the same as in the initial state to within 2%; 

3) the ratio n + /n~ = l.±0.02, which helps constrain the 
isospin asymmetry. 

This last ratio appears redundant, as we already indepen- 
dently use the yields of ir + and ir~ . These yields have 
a large systematic error and do not constrain their ratio 
well, and thus the supplemental constraint is introduced, 
since SHARE allows for the isospin asymmetry effect. 

The successful description of particle yields within the 
SHM obtained for a single chemical freeze-out condition 
produces, as a first result, the model parameters in the 
process of y 2 minimization: the (chemical) freeze-out 
temperature T, the baryon fiB and hyperon /is chemi- 
cal potentials. We obtain and present results at three 
chemical condition alternatives, the chemical equilibrium 
(dashed lines, red online), strange quark non-equilibrium 
with phase space occupancy 7 S ^ 1 (dotted lines, vio- 
let online), and the light quark flavor yield (full) non- 
equilibrium model including "f q =/= 1 (solid lines, blue on- 
line). In our approach, 4 to 6 parameters confront, in a 
systematic fashion, 11 yields and/or ratios and/or con- 
straints containing one redundancy. 

The results we present for the model parameters, in 
Fig. |21 have all above 85% confidence level (we do not 
present the low centrality chemical equilibrium results 
as these do not satisfy this criterion). In table [I] we 
present all these results with precision which should help 
reproduce particle yields when required. In general, the 
reader should not expect to reproduce the last fourth 
digit shown. We also present, in table |]J along with A, 
and in same precision, the volume normalization factor 
dV/dy, which is required to obtain the particle yields. We 
further note that only one of the three models presented 
is applicable, and thus, the variation of parameters be- 
tween these models must be seen as the sensitivity of the 
data to their physical relevance. 

A notable feature, in Fig.[2J is absence of centrality 
features in temperature T and chemical potentials (Ib,Si 
except for the most peripheral centrality, and up to the 
variation which can be associated with fluctuation in the 
data sample and/or determination of the minimum of 
X 2 - The deviation, at the most peripheral centrality bin 
from trends set by other results, could be an indication 
of the change in the reaction mechanism for which we are 
looking. 
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FIG. 2: (color online) From top to bottom: temperature T, 
light quark phase space occupancy 7 9 , the ratio of strange to 
light quark phase space occupancies 7s/7<j and the chemical 
potentials (B for baryochemical fiB and S for strangeness fis) 
as function of centrality (average participant number A). The 
lines connect the results obtained at each bin center A, and 
have not been smoothed in order to show result fluctuations: 
1) Full chemical non-equilibrium model — (blue online) solid 
lines; 2) Strangeness chemical non-equilibrium model — (vi- 
olet online) dotted lines; 3) Chemical equilibrium model — 
(red online) dashed lines. 



Independent of the chemical (non-)equilibrium as- 
sumption, the baryochemical potential fis = 25 ± 1 
MeV across the 10 centrality bins. Similarly, we find 
strangeness chemical potential fj,s = 5.5 ±0.5 MeV 
(related to strange quark chemical potential fi s = 
/zs/3 — us)- The freeze-out temperature is for the semi- 
equilibrium and equilibrium model about 10% greater 
than the full chemical non-equilibrium freeze-out. The 
most notable variation, in Fig-El is the gradual increase 
in strangeness phase space occupancy "fs/jq and thus 
strangeness yield with collision centrality. This effect was 



TABLE I: Fitted statistical parameters, for each central value 
of centrality expressed in terms of participant number A, as 
defined by PHENIX; for the three models in sequence: chem- 
ical non-equilibrium model, chemical semi-equilibrium, and 
chemical equilibrium. For dV/dy the unit is fm 3 , for T, /Mb, A*s 
the unit is MeV, all other quantities are dimensionless. 
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predicted and originates in an increasing lifespan of the 
fireball 0] . The over-saturation of the phase space has 
been also expected due to both, the dynamics of expan- 
sion |17|. and/or reduction in phase space size as a parton 
based matter turns into HG 18] . This latter effect is also 
held responsible for the saturation of light quark phase 
space 7 g — > e m ^/ 2T _ A systematic increase of 7 S with 
collision centrality has been reported for several reaction 
energies |Tg| . 

Using SHARE, we find T = 155±8 MeV for the chem- 
ical equilibrium and strangeness non-equilibrium freeze- 
out. The error is our estimate of the propagation of the 
systematic data error, combined with the fit uncertainty; 
the reader should note that the error comparing central- 
ity to centrality is negligible. This result for T is in mild 
disagreement (1.5 s.d.) with earlier chemical equilibrium 
fits 0,H3|. This, we believe, is due to: 1) differences in 
data sample used, specifically, the hadron resonance pro- 
duction results used provide a very strong constraint for 
the fitted temperature, and 2) the more complete treat- 
ment by SHARE of hadron mass spectrum and of heavy 
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resonances multi-particle decays. 

We find two regimes of hadronization temperature: 

a) for the chemical equilibrium case, and for the chemi- 
cal semi-equilibrium case the hadronization temperature 
T = 155 ± 8 MeV is right at the value of the cross- 
over temperature for lattice QCD with 2+1 dy namical 
flavors results for /ib — * in Refs. [2l| and [23. These 
two groups apply different methods and approximations, 
their results imply that at RHIC conditions the tempera- 
ture of the cross-over from QGP to HG in chemical equi- 
librium is near to T = 160 ± 5 MeV. 

b) The chemical non-equilibrium model requires for in- 
ternal consistency a fast hadronization process, for which 
the sole known mechanism is super cooling of the fireball 
due to rapid transverse expansion j2^| . Supercooling im- 
plies that the hadronization temperature is below phase 
cross-over boundary, and the value T = 141 ± 7 MeV, we 
have found in the full chemical non-equilibrium fit, is a 
result internally consistent with this reaction mechanism. 
These issues are discussed further in the context of the 
search of phase threshold as function of reaction energy, 
seeRef. [l3- 

For the chemical non-equilibrium sudden hadroniza- 
tion reaction picture, we further expect that the 
hadronization geometry is highly non-homogeneous, cor- 
responding to, e. g., fingering breakup of the fireball. This 
makes it possible that the thermal and chemical freeze- 
out condition coincide 0> HH , as the surface to volume 
area is large, emitted particles having a small probabil- 
ity to rescatter. For this reason, as noted before, there 
is no need to discuss changes in observable hadron reso- 
nance populations which may arise in post-hadronization 
scattering processes. 



III. PHYSICAL PROPERTIES OF THE 
FIREBALL AT HADRONIZATION 

Given the statistical parameters, we can evaluate the 
yields of particles not yet measured and obtain the rapid- 
ity yields of entropy, net baryon number, net strangeness, 
and thermal energy as function of centrality, shown in 
Fig.|3| We find, in the most central reaction bin, 14.9±1.5 
baryons per unit rapidity interval, a rather large baryon 
stopping in the central rapidity domain. 

The entropy yield is reaching dS/dy — 5000 ± 500. 
This compares to the estimate made recently by Pal and 
Pratt H3 who find dS/dy = 4451 ± 445 for the most 
central 130 GeV reactions. One should note that the 
smoothness of the results, presented in Fig.|3|as function 
of A, is directly related to the smoothness of experimen- 
tal data, which determines in some cases nearly directly 
these observables. For example, the hadron multiplicity 
per rapidity unit is directly related to the entropy yield. 
Thus, the experimental yield error of 10% is directly the 
error of entropy dS/dy. If instead, we constructed this 
error from partial errors in the statistical parameters, the 
implicit cancellations would be hard to realize. For this 
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FIG. 3: (color online) From top to bottom: entropy, net 
baryon, strangeness and thermal energy yield per unit of ra- 
pidity, as a function of centrality. Lines (nearly overlapping) 
are coded as in figure |21 for the three chemical models. 



reason, we do not state these errors in Fig. [3 all these 
micro canonical quantities are directly relate to particle 
yields and their error is at 10% level. 

Except for the chemical equilibrium model (which has 
a restricted centrality range of validity, for it fails to fit 
the peripheral collisions), the rise of strangeness yield 
with centrality is faster than the rise of baryon number 
yield: (ds / dy) / (dB / dy) = s/B is seen in the top panel 
of Fig.0J For the most central reactions, we reach s/B = 
9.6 ± 1. We also note fluctuations in the trend of the 
results of the magnitude expected from the experimental 
data error. 

The increase, with A, of per baryon specific strangeness 
yield indicates presence of a production mechanism act- 
ing beyond the first collision dynamics. This new mecha- 
nism must benefit from the increased size, or more appro- 
priately, increased life span of the larger reaction system. 
We have shown earlier that the thermal gluon fusion to 
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FIG. 4: (color online) From top to bottom : strangeness 
per net baryon s/B and strangeness per entropy s/S, as a 
function of centrality. Lines are coded as in figure [5] for the 
three chemical models. 



strangeness will have just this behavior [l6| . 

While entropy production occurs predominantly dur- 
ing the initial parton thermalization phase, thermal 
strangeness production requires presence of thermal, mo- 
bile gluons, being driven by thermal gluon fusion [2fj| . 
Thus, strangeness production follows in time the entropy 
production, and is strongest at the highest available tem- 
perature considering the strangeness mass threshold. 

To understand this better, we show, in the bottom 
panel of Fig.0] strangeness per entropy s/S, as function 
of centrality. We hnd a smooth transition from a flat 
peripheral behavior where s/S ^ 0.02, to a smoothly in- 
creasing s/S reaching s/S ~ 0.03 for most central reac- 
tions. This rise, occurring for A > 20, indicates the onset 
of an additional strangeness production mechanism. 

As the system size increases, the time span, during 
which thermal strangeness production is effective, is in- 
creasing allowing strangeness to approach chemical equi- 
libration in the parton phase. A back of envelope es- 
timate of the ratio of s/S, in a chemically equilibrated 
parton plasma, yields the value seen in Fig.|3j for large 
A 0. We thus understand the rise and magnitude of 
s/S ratio in terms of the expected reaction mechanisms 
in the deconfined phase, and the chemical saturation of 



this ratio. Chemical saturation in QGP, i.e., 7 
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1, 



implies, due to higher strangeness phase space content in 
QGP than in hadron matter, that 7^ > 1. 

We evaluate now the bulk properties of the fireball 
shown in Fig.[S] To obtain the pressure P, the energy 
density e = (dE/dy) /(dV/dy) and the entropy density 
a = (dS/dy) /(dV/dy), we sum, using the SHARE data 
base for hadron resonances, all partial particle contri- 
butions using relativistic expression, see 10]. Since we 
fitted six particle rapidity yields, the global fitted yield 
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FIG. 5: (color online) From top to bottom: pressure P, 
energy density e = E/V, entropy density S/V and E/TS, as 
a function of centrality. Lines are coded as in figure [5] for the 
three chemical models. 



normalization factor dV(A)/dy is reliable, up to the sys- 
tematic error of about 10% inherent in the PHENIX ra- 
pidity yield data. 

When we consider ratios of two bulk properties, e.g., 
E/TS, the data and fit fluctuations cancel out, and the 
results are in general much smoother. Moreover, such 
ratios do not depend on the identification in the nor- 
malization factor dV/dy with the statistical model vol- 
ume. For example, when one considers hadrons of finite 
proper volume |28j . there is a correction factor. Hence, 
the smooth and precise values for E/TS which are not 
expected to change in their interpretation offer a chal- 
lenge for structure model of hadronization, and will need 
to be addressed quantitatively. We note, as an example, 
that a quark matter system consisting of thermal mass 
quarks with m = aT, a ~ 2 for non-equilibrium and 
a = 4 for equilibrium, will yield just the result seen in 
the bottom panel of Fig.0 

We note, in Fig.0 that the chemical non-equilibrium 
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system is much denser at hadronization |29j . The entropy 
density, assuming chemical non-equilibrium, is by more 
than a factor 2 larger compared to the chemical equilib- 
rium model. It is for this reason that the ratio E/TS, in 
Fig-El is f° r t ne non-equilibrium case significantly smaller 
than it is when assuming equilibrium. 

IV. FINAL REMARKS 

We have presented a comprehensive analysis of soft 
hadron yields at ^/snn = 200 GeV as function of central- 
ity We have obtained the statistical hadronization pa- 
rameters that describe the data, and evaluated the phys- 
ical properties of the hot fireball at hadronization. Our 
analysis included, aside of 'stable' hadron PHENIX data 
(tt ± , K ± ,p,p), also the STAR K* and <f> yields. For the 
latter, we also presented a comparison analysis in order 
to resolve a discrepancy between STAR and PHENIX 
0-yield result. We found that, for A £ 20, the statisti- 
cal parameters do not vary with centrality, with the ex- 
ception of strangeness quark occupancy, 7 S . Already for 
A £ 20, the system properties are approaching those seen 
for the greatest available A ~ 350. However, for A < 20 
there is a significant change in the physical and statisti- 
cal properties of the fireball. The chemical equilibrium 
description of the experimental results here considered is 
not possible for these small fireballs. 

We are aware of two prior efforts to explore statistical 
parameters, but not bulk properties, as function of colli- 
sion centrality. At ^/snn = 130 GeV [3(|, the analysis re- 
mains inconclusive in view of the limitations of the exper- 
imental data. At ^/snn = 200 GeV, an analysis assuming 
chemical semi-equilibrium |3l| shows trends comparable 
with those we found. However, our hadronization tem- 
perature T and phase space occupancy 7 S are anchored 
by the K*(892) and (j> yields, and hence, we obtain for 
the chemical semi-equilibrium a smaller value of T, and 
therefore a greater 7„ which for most central reactions 
clearly exceeds unity. 

In this work, we have obtained the centrality depen- 
dence of the hadronization pressure, entropy density and 



thermal energy density, which are of the magnitude ex- 
pected: for the chemical non-equilibrium model the pres- 
sure saturates at P = 92 MeV/fm 3 , a typical value we 
are used to from the bag model of hadrons to be the 
vacuum pressure. The energy density of 0.5 GeV/fm 3 
is in accord with lattice results for the energy density 
of matter subject to the phase transformation to/from 
deconfincment . 

We have shown that the fireball strangeness content is 
increasing fastest with increasing centrality, beating out 
in the competition both, the stopped net baryon number, 
and the produced entropy. We have shown that the ratio 
s/S saturates near the value of chemically equilibrated 
QGP phase. 

In our opinion, the most remarkable finding of this 
study is the recognition that the statistical parameters, 
and thus also the bulk properties of dense matter fireball 
created at RHIC, do not depend on the size of the system 
for A ^ 20, where A is the number of reaction partici- 
pants. There is a rapid adjustment in the value of statis- 
tical parameters fitted: temperature is dropping and ~/ q 
is increasing, which changes combine to yield a doubling 
in the density of the hadronizing fireball for A ^ 20. In- 
terpreting the high density phase at hadronization as the 
deconfined state, our results can be interpreted to indi- 
cate that at RHIC energy scale the quark liquid phase is 
formed for A £ 20. 
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